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Abstract 

We consider the existence of invariant manifolds in phase space governing reaction dynamics in 
situations where there are no saddle points on the potential energy surface in the relevant regions 
of configuration space. We point out that such situations occur in a number of important classes 
of chemical reactions, and we illustrate this concretely by considering a model for transition state 
switching in an ion-molecule association reaction due to Chesnavich (J. Chem. Phys. 84, 2615 
(1986)). For this model we show that, in the region of configuration space relevant to the reaction, 
there are no saddle points on the potential energy surface, but that in phase space there is a 
normally hyperbolic invariant manifold (NHIM) bounding a dividing surface having the property 
that the reactive flux through this dividing surface is a minimum. We then describe two methods 
for finding NHIMs and their associated phase space structures in systems with more than two 
degrees-of-freedom. These methods do not rely on the existence of saddle points, or any other 
particular feature, of the potential energy surface. 
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I. INTRODUCTION 



Critical points of the potential energy surface have played, and continue to play, a sig- 
nificant role in how one thinks about transformations of physical systems^. The term 
'transformation' may refer to chemical reactions such as isomerizations^— or the analogue 
of phase transitions for finite size systems^^. A comprehensive description of this so- 
called 'energy landscape paradigm' is given in ref. O. The energy landscape approach is an 
attempt to understand dynamics in the context of the geometrical features of the potential 
energy surface, i.e., a configuration space approach. However, the arena for dynamics is 
phase spaced— , and numerous studies of nonlinear dynamical systems have taught us that 
the rich variety of dynamical behavior possible in nonlinear systems cannot be inferred from 
geometrical properties of the potential energy surface alone. (An instructive example is the 
fact that the well-studied and nonintegrable Henon-Heiles potential can be obtained by se- 
ries expansion of the completely integrable Toda system^.) Nevertheless, the configuration 
space based landscape paradigm is physically very compelling, and there has been a great 
deal of work over the past ten years describing phase space signatures of index one saddles^- 
of the potential energy surface that are relevant to reaction dynamics (see, for example, refs 



22 



2J). More recently, index two^— and higher index— saddles have been studied. 



The work on index one saddles has shown that, in phase space, the role of the saddle point 
is played by an invariant manifold of saddle stability type, a so-called normally hyperbolic 
invariant manifold or NHIM— >22. The NHIM proves to be the anchor for the construction 
of dividing surfaces that have the properties of no (local) recrossing of trajectories and 
minimal (directional) flux—. There is an even richer variety of phase space structures and 
invariant manifolds associated with index two saddles of the potential energy surface, and 
their implications for reaction dynamics are currently under investigation^. Fundamental 
theorems assure the existence of these phase space structures and invariant manifolds for 
a range of energy above that of the saddle^. However, the precise extent of this range, as 
well as the nature and consequences of any bifurcations of the phase space structures and 
invariant manifolds that might occur as energy is increased, is not known and is a topic of 
current investigation^. 

While work relating phase space structures and invariant manifolds to saddle points on 
the potential energy surface has provided new insights and techniques for studying reaction 
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dynamics^ — } it certainly does not exhaust all of the rich possibilities of dynamical phe- 
nomena associated with reactions. In fact, recent work has called into question the utility 
of concepts such as the reaction path and/or transition state^r— . Of particular interest for 
the present work is the recognition that there are important classes of chemical reaction, 
such as ion-molecule reactions and association reactions in barrierless systems, for which the 
transition state is not necessarily directly associated with the presence of a saddle point on 
the potential energy surface (or even the amended potential, which includes centrifugal con- 
tributions to the energy^ 1 ^). The phenomenon of transition state switching in ion- molecule 
reactions^— provides a good example of the dynamical complexity possible in such systems. 

The lack of an appropriate critical point on the potential energy surface with which to 
associate a dividing surface separating reactants from products in such systems does not 
however mean that there are no relevant geometric structures and invariant manifolds in 
phase space. In this paper we discuss the existence of NHIMs, along with their stable 
and unstable manifolds and associated dividing surfaces, in regions of phase space that do 
not correspond to saddle points of the potential energy surface. After presenting a simple 
example motivated by Chesnavich's model for transition state switching in an ion-molecule 
association reaction^, we describe a theoretical framework for describing and computing 
such NHIMs. Like the methods associated with index one and two saddles, the method we 
develop for realizing the existence of NHIMs is based on normal form theory; however, rather 
than normal form theory for saddle-type equilibrium points of Hamilton's equations (which 
are the phase space manifestation of index one and two saddles of the potential energy 
surface), we use normal form theory for certain hyperbolic invariant tori. The hyperbolic 
invariant tori (and their stable and unstable manifolds) alone are not adequate, in terms 
of their dimension, for constructing NHIMs that have codimension one stable and unstable 
manifolds (in a fixed energy surface). However, by analogy with the use of index one saddles 
to infer the existence of NHIMs (together with their stable and unstable manifolds, and other 
dividing surfaces having appropriate dimensions), these particular hyperbolic invariant tori 
can likewise be used to infer the existence of phase space structures that are appropriate for 
describing reaction dynamics in situations where there is no critical point of the potential 
energy surface in the relevant region of configuration space. 

Section |TT] discusses our simplified version of Chesnavich's model for transition state 
switching^. For this 2 DoF system, we exhibit a NHIM (in this case, an unstable peri- 
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odic orbit) that is the rigorous dynamical manifestation of the mininimal flux surface of 
variational transition state theory^— . In Section III II we describe a (time-dependent) nor- 
mal form based approach for finding such NHIMs in phase space. In particular, we present 
two variations of the method. In Section IHIAI we consider systems where (to leading order) 
the system can be separated into a two degree-of-freedom (DoF) subsystem and a collection 
of decoupled "bath modes" . We assume that there exists a hyperbolic periodic orbit in the 
2 DoF subsystem, and show that this can be used to construct a hyperbolic torus for the 
full system. We then show that this hyperbolic torus implies the existence of a NHIM, with 
stable and unstable manifolds that are codimension one in the energy surface. Appropriate 
dividing surfaces can then be constructed using the NHIM and the normal form Hamilto- 
nian. In Section IIII Bl we describe a method which requires knowledge of the appropriate 
hyperbolic invariant torus from the start. The advantage of the first method is that it is 
more intuitive and can exploit the considerable number of methods for locating hyperbolic 
periodic orbits in 2 DoF Hamiltonian systems. Method two is more general, but at present 
there are few techniques available for locating hyperbolic invariant tori of the appropriate 
dimension in general N DoF Hamiltonian systems that are not perturbations of integrable 
systems. Section [TV] concludes. 
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II. A MOTIVATING EXAMPLE: VARIATIONAL TRANSITION STATE FOR A 
MODEL BARRIERLESS REACTION 

A. Introduction 

The conventional approach to variational transition state theory (VTST) for barrierless 
reaction proceeds by minimizing the reactive flux with respect to variation of some reaction 
coordinate chosen a priori^—. The value of the reaction coordinate so determined is there- 
fore the location of a flux bottleneck, which is identified with the transition state for the 
particular association reaction. An invariant phase space characterization of such variation- 
ally determined dividing surfaces is highly desirable; for N = 2 DoF, such transition states 
will presumably be associated with unstable periodic orbit dividing surfaces (PODS)^r— , 
or, more generally, with NHIMs (N > 2 DoF)^o. 



B. Model Hamiltonian 



We consider a highly simplified model for a barrierless association reaction (cf. ref. l43f ). 
The system has 2 DoF: a radial coordinate r, identified as the reaction coordinate, and a 
coordinate s describing vibrations transverse to the reaction coordinate. The radial poten- 
tial has the character of a long-range attractive ion-neutral interaction, while the potential 
transverse to the reaction coordinate is harmonic. The system is nonseparable by virtue of 
a dependence of the harmonic oscillation frequency on the coordinate r: 



2 2 r 4 2 v 

We take u)(r) to have the form 



u(r) = u e-P r (2) 

so that, for f3 > 0, the transverse vibration stiffens as r decreases. 
A contour plot of the potential 

v(r,s) = -— + -w(r)V (3) 
r 4 2 

for parameter values a = l,/3 = l, Wo = 8is shown in Fig. [TJ 
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C. Locating the bottleneck 



For the 2 DoF Hamiltonian ([T]), we can compute the action of the transverse vibrational 
mode as a function of the coordinate r at fixed energy E: 



As r decreases, there are two competing tendencies: 

• Decreasing r increases the amount of energy in the oscillator degree of freedom, thereby 
tending to increase the action. 

• Decreasing r increases the frequency oj(r), tending to decrease the action of the tran- 
verse mode. 

The competition between these two trends can therefore result in the existence of a 
minimum in the action as a function of r: see Figure EJ 

The minimum of the action as a function of r corresponds to an extremum of the sum of 
states (phase space area) or flux as a function of r, and hence is interpreted as a bottleneck. 
In the variational transition state approach, the transition state for association is then 
located at the value r = r* corresponding to minimum flux. This bottleneck corresponds to 
an inner or "tight" transition state^ 1 ^. 

D. Intrinsic characterization of variational TS: PODS 

The formulation of VTST outlined above for the model association reaction is unsat- 
isfactory in that the minimum flux bottleneck so determined has no intrinsic dynamical 
significance. It is natural to seek a dynamical, phase space based characterization of the 
variational TS. For 2 DoF systems, transition states are identified as PODS^— . The invari- 
ant object defining the TS is a hyperbolic (unstable) periodic orbit; the ID periodic orbit 
forms the boundary of a 2D dividing surface on the 3D energy shell in phase space, which is 
the transition state. The minimal flux (local no-recrossing) property of the TS follows from 
the principle of stationary action^. 

A search using the model potential ([3]) reveals the presence of a PODS in the vicinity 
of r ~ r* (the periodic orbit at E — 1 is shown in Fig. [T])£i. This PODS is the rigorous 
dynamical realization of the variational TS in this simple case. 




(4) 
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III. LOCATING NHIMS WHEN THERE ARE NO (RELEVANT) SADDLES ON 
THE POTENTIAL ENERGY SURFACE 



In this section we describe two methods for locating NHIMs of the type discussed in 
the previous section. These methods are inherently phase space approaches, based on the 
existence of a hyperbolic invariant torus solution of Hamilton's equations. Normal form 
theory for hyperbolic invariant tori can be used to provide "good coordinates" for computing 
explicit formulae for a NHIM, its stable and unstable manifolds, and dividing surfaces in the 
phase space vicinity of the hyperbolic invariant torus on which we base our method, in much 
the same way that it is used to compute similar objects associated with index one saddles 
of the potential energy surface^ 1 ^ 1 ^. A large literature for normal form theory associated 
with invariant tori of Hamilton's equations has been developed over the past twenty years, 



see ref. 



53 



for an overview, and ref. 



54l for a survey of the issues associated with bifurcation 



of tori in Hamiltonian systems. For our present purposes we use the results contained in 



ref. 



55 



which explicitly discusses the relevant normal form and also clarifies the issue of 
"hyperbolicity" of tori in Hamiltonian systems (concerning which there had previously been 
some confusion in the literature). 



A. Method 1: a relevant 2 DoF subsystem can be identified at leading order 

Consider a Hamiltonian of the following form: 

H = y + y + V{r,s) + - J^CJi {uj + vf) + /(r,S,Mi, . . . ,U n _ 2 ,Pr,Ps,Vi, . . . , V n _ 2 ) (5) 

i=l 

where f(r, s,Ui,..., u n _2,Pr,Ps> vi, ■ ■ ■ , Vn-2) is a t least order 3, denoted 
O s (r, s, u\, . . . , u n -2iPr,Psi v ii ■ ■ ■ j v n-2)- hi general this term serves to couple all of 
the variables, but we have written the Hamiltonian in such a way that we can identify a 
clearly defined 2 DoF subsystem, on which we make the following assumption: 

Assumption: The 2 DoF subsystem defined by the Hamiltonian: 

n = ^ + ^ + V(r,s), (6) 

has a hyperbolic periodic orbit, denoted V = (r(t), s(t),p r (t),p s (t)). 
To construct a NHIM for (jSJ) we proceed as follows. 

7 



Step 1: Transform the 2 DoF subsystem to normal form in a neighborhood of 
the periodic orbit. 

we can find an invertible transformation T (as smooth as the Hamil- 



Following ref. 
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tonian) defined in a neighborhood of the periodic orbit 

T : (r,s,p r ,p s ) ^ T (r, s,p r ,p s ) = (l,0,x,y) (7a) 
Tq" 1 : (l,0,x,y) i-» T -1 (J, 0, x, y) = (r,s,p r ,p s ) (7b) 

such that the 2 DoF Hamiltonian takes the form: 

K = cul + Xxy + 2 (I) + 3 (I, x, y), (8) 

where we can take u, A > 0. Of course, nontrivial calculations are required in going from 
(jSJ) to dHJ). In particular, after the hyperbolic periodic orbit is located, a time-dependent 
translation to "center" the coordinate system on the periodic orbit must be carried out; the 
resulting Hamiltonian is then Taylor expanded about the origin (i.e., the periodic orbit), 
a Floquet-type transformation constructed to make (to leading order) the dynamics in the 
normal direction to the periodic orbit constant (i.e., A is constant in (jSJ)), then, finally, 
Hamiltonian normal form theory is applied to the result. Details of the methodology for 



carrying out this procedure for specific examples are described in refs l56ll57l . For the purposes 
of demonstrating the existence of a NHIM, we only need to know that such transformations 
can in principle be carried out. 

Step 2: Use the normal form transformation for the 2 DoF subsystem to rewrite 



We use the normal form transformation of the 2 DoF subsystem to express (jSJ) as follows: 

n-2 



n-2 

H = ujI + Xxy + - oji (uf + vf) 



i=l 



+ F(u u u n _ 2 , vi,..., v n - 2 , 1, 0, x, y) + 2 CO + 3 (I, x, y) (9) 

where 

F(ui, u n -2i vx,..., v n - 2 , 1, 9, x, y) = 

f(r(I, 9, x, y), s(I, 9, x,y),u h ..., u n _ 2 ,Pr(I, 0, x, y),p,(I, 9, x, y), v x , ■ ■ - , v n _ 2 )- (10) 

Step 3: Use Hamiltonian (jUj) to conclude the existence of a NHIM. 
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Let 



I=h v ? + #) i O = tan^f-). (11) 



We then rewrite ([9]) as 
K = Xxy + |(w 2 + z 2 ) + - ^ Ui (u 2 + W 2 ) 

i=l 

+ F(u u ii n _ 2 , «i, • • • , v„_2, /(tu, z),0{w, z),x, y) + 2 {I(w, z)) + 3 (J(iw, z), x, y). (12) 



By construction, "H(O) = 0. Neglecting higher order terms in ffT2"l) . we obtain: 

n-2 



^ n-2 

^trunc = + ^(U> 2 + Z 2 ) + - ^ CJ; [u] + W 2 ) . (13) 



j=l 



If we set x = y = 0, then on the energy surface "Htmnc = /i > 0, 

|(io a + ^) + i£ Wi W + t,?)=fc>0 ) (14) 
l=\ 

is a normally hyperbolic invariant 2n — 3 sphere in the 2n-dimensional space with coordinates 
(ui, . . . , u n -2, V\, . . . , v n -2, w, z, x, y), having 2n — 2 dimensional stable and unstable mani- 
folds in the In — 1 dimensional energy surface. The persistence theory for NHIMs implies 
that this manifold persists when the higher order terms are added for energies sufficiently 



close to h — (cf. ref. l30f ). 

At this point we are in a position where normal form theory can be used on ( TT2|) to 
construct a new set of coordinates (the normal form coordinates) in which ( fl~2l) assumes 
a particularly simple form that results in explicit formulae for the NHIM, its stable and 
unstable manifolds, and dividing surfaces between regions of the phase space corresponding 
to reactants and products. The normal form algorithm also provides the transformation from 
the original physical coordinates and its inverse, and this allows us to map these surfaces 



back into the original physical coordinates, as described in ref. 



24 



B. Method 2: A hyperbolic torus of dimension n — 1 can be located in an n DoF 
system. 

The advantage of method 1 is that it can make use of extensive prior work on locat- 
ing periodic orbits, and determining their stability, in 2 DoF systems. In method 2 the 
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starting point is knowledge of the existence of a hyperbolic torus of dimension n — 1 in 
the 2n-dimensional phase space, with the frequencies on the torus satisfying a diophantine 



condition^. In this case, it follows from ref. 



55 



that an invertible transformation of coordi- 



nates, valid in a neighborhood of the torus, can be found where the system has the following 
form: 

H = uili + . . . + oj n -iI n -i + Xxy + 02(A) • ■ • , In-i) + . . . 

y)- (is) 

Now let 

k = \[u\ + «?), 0i = tan" 1 f^J , i = 1, . . . , n - 1. (16) 
In terms of these coordinates the Hamiltonian (1171) has the form: 



n-1 



H = -^2u i (u* + Vi) + \xy + 2 (u 2 1 + vl, . . . , u^_i + v^i) + 3 («i + «i, • • • » «n-i+ v S-n ^2/) 

i=l 

(17) 

We now proceed exactly as for method 1. Neglecting the terms in (ITT]) of order 3 and 
higher gives: 

- n— 1 

-^Wi(M? + T;?) + Axj/. (18) 



n-1 

Ft, 



[ trunc 2 

i=l 

If we set x = y = 0, then on the energy surface H trunc = h > 0, 

^ n— 1 

-^^K 2 + ^ 2 ) = / i >0, (19) 

i=l 

is a normally hyperbolic invariant 2n — 3 sphere in the 2n-dimensional space with coordinates 
(ui, . . . , u n -\, V\, . . . , v n -i, x, y) having In — 2 dimensional stable and unstable manifolds in 
the 2n — 1 dimensional energy surface. The persistence theory for NHIMs implies that this 
manifold persists when the higher order terms are added for energies sufficiently close to 
h = 0. 

Method 2 relies on first finding an appropriate hyperbolic invariant torus, and in general 
the existence of this object will need to be verified numerically. There has been a great deal 
of activity developing such numerical methods in recent years. See, for example, refs 



59 



66 



Finally, it is also worth noting the difference between the h = limit for index 1 saddles 
on the potential energy surface compared with the torus case in phase space. In the former 
case the NHIM shrinks down to a point on the potential energy surface (i.e. configuration 
space), while in the latter case it shrinks down to a torus in phase space. 
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IV. SUMMARY AND CONCLUSIONS 



In this paper we have exhibited a normally hyperbolic invariant manifold (NHIM^, in this 
case a PODS^r— , or unstable periodic orbit) defining a flux bottleneck in a simple model 
of an ion-molecule reaction, for which there is no associated critical point of the potential 
energy surface. We have also developed a theoretical framework for showing the existence 
of such NHIMs. Two methods were described that are in principle suitable for computing 
such phase space objects in the multidimensional case. 

For index one saddles a software package has been developed that allows one to com- 
pute the normal form associated with the corresponding saddle-center-. . .-center stability 
type equilibrium point to high order for multi-dimensional systems, with control over the 
accuracy^. Accuracy is assessed by a battery of tests, and specifying an accuracy may 
affect the order of the normal form that can be computed as well as the dimensionality of 
the system that can be treated. In general, these issues must be analyzed on a case-by-case 
basis. Nevertheless, the normal form, and most importantly the transformation from the 
original physical coordinates and its inverse, allow us to realize the NHIM, its stable and 
unstable manifolds, and dividing surfaces between regions of the phase space corresponding 
to reactants and products. Moreover, flux through the dividing surfaces can be computed 
as an integral over the NHIM, and the normal form coordinates provide a natural way of 
selecting distributions of initial conditions of trajectories on the dividing surfaces to com- 
pute gap times^^. It should be possible to develop similar software for computing normal 
forms associated with hyperbolic tori of the type discussed above, and one expects that the 
normal form will, similarly, allow one to realize phase space structures relevant to reaction 
dynamics as well as to compute fluxes and sample distributions of initial conditions. 

This program has yet to be carried out, but the essential computational elements of the 



approach can be found in refs 



56 



571 . Such capabilities would be very useful for the study of 



reaction dynamics in multimode systems exhibiting transition state switching, for example 
ref. 42. 
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Figure captions 



FIG. 1: Contour plot of the model potential of eq. ([3]). Parameter values a — 1, /3 — 1, 
ujq = 8. Superimposed on the potential contours is a PODS defining a "tight" transition 
state, computed at energy E — 1. 

FIG. 2: Action I(E,r) (eq. (j4])) as a function of r, computed at energy E = 1. 
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FIGURE 2 
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